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Abstract: A method for clean determination of the unitarity 
angles a and 7 is proposed that uses only direct CP violation 
^ \ and does not require any time dependent measurements. The 

method takes advantage of helicity amplitudes for B u , Bd and B s 
^ ■ decay to two vector mesons and can be used, at any B-facility, 

O . in conjunction with a large number of modes. It also allows for 



experimental tests of theoretical approximations involved. 



Considerable progress has recently been made in experimental determina- 
tion of the angle (3 of the unitarity triangle and improved measurements are 
D \ expected in the near future [|], |2|, However, even with increased precision 

^ ; in the extraction of f3, the ability to test the CKM ||J description of CP vio- 

lation in the Standard Model (SM) will be limited as the existing tests rely 
^ ■ on theoretical calculation of hadronic matrix elements which still have con- 

siderable uncertainties ||. Therefore clean determinations of all three angles 
(a, (3, 7) of the unitarity triangle is important to facilitate precision tests of 
the SM and search for new CP-odd phase(s) due to physics beyond the SM. 
Currently the asymmetric B-factories are performing remarkably well and in 
the future hadronic B-machines may produce 1-3 orders of magnitude more 
B-mesons. Thus, it is very important to devise methods for clean extraction 
of the angles that can be used at all B-facilities. Furthermore, since time de- 
pendent measurements involving B s mesons are extremely challenging due to 
the high frequency of B s — B s oscillations, it would also be helpful if methods 
could be devised that do not require such time dependent information. 
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Motivated by these considerations we propose a method with the novelty 
that it uses only direct CP-violation and does not involve any use of time 
dependent measurements. Potentially all types of B-mesons (B u , B d , B s ) can 
be used at all B-facilities. The method also allows for experimental tests of 
the key theoretical approximations that are involved. Specifically, we propose 
determination of the angles a and 7 through a study of the interference of 
tree and penguin amplitudes in decays of B-mesons to two vector particles. 

There have been several methods proposed to extract these angles. For a 
one can consider oscillation effects in B° — > 7r + 7r~ although one must account 
for the penguin through isospin analysis || by observing B° — > 7r°7r°. Since 
the branching ratio to 7c°tt° is expected to be small and hard to observe, 
it may be preferable to consider B° — > 7r + 7r~7r° where one can also take 
advantage of resonance effects in the Dalitz plot 0; however there may be 
problems in precise modeling of the resonance structure. Another method 
for extracting a from the interference of w-penguins with t-penguins in B° — > 
may overcome the disadvantages of the 2ir and 3ir final states |§ 
although the analogous B s decays are required for the analysis. All three 
of these methods use time dependent CP violation measurements of B°(B°) 
decays driven by the mixing in that system while our method relies on direct 
(time independent) CP violation only. 

The angle 7 may be extracted at the S-factories through the interference 
of b — > ucs and h — > cus || . In the method we discuss here, we will obtain 7 
through the interference of the b — > s penguin and the b — > suu tree and a 
from the interference of the b — > d penguin and the b — > duu tree. 

Our method requires measurements made in two separate decay modes. 
One mode of the form B — > V1V2, together with B — > V1V2, (B = B u , Bd or 
B s ) which receives both tree and penguin contributions and another mode 
B — > V3V4 which is a pure penguin. The procedure uses helicity amplitudes of 



B — > V1V2, which can be inferred from decay distributions of the vectors [10 



|TT| , O O, . By analyzing this data together with the helicity amplitudes 
for B — > V3V4 one can extract the tree-penguin weak phase difference. From 
a b — * d penguin transition contributing to B — > V1V2, we deduce |^sina| 
from which | sina| may be obtained once |Vtd/T4 s | is known [15[]; from a b — > s 
penguin we can obtain | sin 7]. A list of such modes is shown in Table [l]. Thus 
to implement this method to obtain a one needs to study one mode from the 
a-mode column and one mode from the "pure penguin" column. Likewise 
to obtain 7 one 7-mode and one pure penguin mode is required. Using more 
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than one mode from each column will, of course, increase the analyzing power 
of the method. 

Let us first discuss some of the tree penguin interference modes that are 
suitable in our method for extracting a. In charged B decays, the relevant 
mode is — > p ± u. In this case the tree graph is color allowed which we 
find is advantageous in terms of statistical power. Since the charge of the p 
indicates the flavor of the initial state, no other tagging is required. However, 
this mode has two 7r°, in the final state originating from decays of p^ and uj 
perhaps making it difficult for hadronic machines. 

From Table [I] we see see that for an analogous extraction of a via decays 
of B° we may use any one of the modes B° — > p°p°, p + p~, u;p° or ujuj. 
These modes require tagging at production and C-odd interference terms are 
degraded somewhat by time integration. Also the tree is color suppressed; 
however, in the case of p°p° it leads to a final state that does not contain 
any 7r° and so may be easier to implement at hadronic B machines. The 
mode B° — > p + p~ has the advantage that any electro-weak penguin (EWP) 
contribution is color suppressed. 

Decays of B s usable for a are B s — > K*°p° or B s —> K* + p~ . Both are self 
tagging, though in the first case only if K*° — > K + n~ . The first mode also 
has the advantage that the final state contains no tt°. The final state K* + p~ , 
while suffering from at least one 7r° in the final state, has the advantage that 
possible contamination from the EWP is color suppressed. 

Likewise, our method can also be used to extract 7 through the use of 
b — > s penguins. In B° the candidate modes are B° — > K* + p~ , B° — > K*°p° 
and B° — > K*°u; note that B° — > K*°p° has the advantage of no 7r° in the 
final state. For charged B's again there are two such modes: B + — > K* + p° 
and B + — > K* + u, either one of which could be used to obtain 7; B + — > K* + p° 
is free of 7r°. 

Our method also needs input from one pure penguin mode; this should be 
a, b —>■ s penguin since we need the component with an intermediate u quark 
to be small. For pure penguin modes we can use: B + — > <pK* + , B° —>■ <pK*°, 
B+ _> K*°p+ 1 B s -> 00 or B s -> K*°K*°. Likewise the decays B° -> p°K*° 
and B° — > uK*° can also be used since the tree contribution is color and 
CKM suppressed. In Table |l] as well as listing all the modes mentioned, we 
indicate which have color suppressed EWP (underlined); which have color 
allowed tree contributions (parentheses) and which have 7r°-free final states 
[square brackets]. 
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Table 1: The relevant B — > VV modes originating from each kind of B 
meson are shown. In the "a-mode" column are tree-penguin interference 
modes which are sensitive to a; likewise the "7-modes" are sensitive to 7. 
The "pure penguin" modes proceed through b — ► s penguin processes only. 
The underlined modes have a color suppressed EWP contribution; the modes 
enclosed in parentheses have color allowed tree contributions while the modes 
enclosed in square brackets have 7r° free final states. 



5-meson 


a-mode 


7- mode 


Pure Penguin 


B + 






[<j)K* + ], K*°p + 


B° 


{p + P~l [p°P% 

UJUJ, p°oj 


(K*+p~) 


[K*°<j>] 


B s 


[K*°p% (K*+p-), 
K*°uj 


see note [16 


[66\. [K*°K*°] 



In our analysis we will use the approximations that: (1) SU(3) is a valid 
symmetry for penguin processes, (2) the effects of the EWP are small, (3) the 
qq pair which arises in a strong penguin does not form a single vector meson 
of the final state. Recall that the EWP are assumed to be small for other 
proposed methods || [7], || as well for extracting a. Note also that for some 
of the modes in Table [l] the EWP is color suppressed. Later we will discuss 
ways to test each of these three approximations. 

To illustrate our method, we now focus on B ± — > p ± u (i.e. V1V2 = p + u) 
where the discussion easily generalizes to the other modes mentioned above. 
The SM amplitude A for this process can be written: 

A = Tv u + P u v u + P c v c + P t v t = fv u + Pv t (1) 

where T is the tree contribution to the amplitude, Pi is the penguin contri- 
bution due to the diagram with an internal quark of type U{ and Vi = V* b Vid- 
Unitarity of the CKM matrix allows us to express the amplitudes in terms of 
T = T + P u — P c which we will refer to as the corrected tree and P = P t — P c , 
the corrected penguin. 
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For this mode it is useful to follow a (non-standard) convention for the 
weak phase where the corrected tree is zero. The amplitudes for B + — > p + uo 
may then be written as: 



A = a + he ia where a = f\v u \ and b = P\v t \ (2) 

and each of the bold face terms represent a set of helicity amplitudes. 

The amplitude for the charge conjugate process in this convention is A = 

a + Uhe~ ia where II indicates parity. 

Even knowing the full amplitudes in these decays we cannot hope to 
determine a. This is evident writing the amplitude: 



A = c + ibsina, A = c — ilTbsina =^> b sin a = (A — IIA)/(2z) (3) 

where c = a + bcosa. While we may extract c and bsina, we clearly can 
not separately obtain sin a without more information [17|. To solve for a 



then, we need two more pieces of information: (1) A means to fix the relative 
weak phase between A and A. (2) A separate normalization of the penguin 
contribution b. 

The measurement of pure penguin amplitudes will provide these remain- 
ing two inputs. The normalization of such a mode will clearly allow the 
determination of sin 2 a while the ratios and phase differences between the 
components of b can be matched with those determined from eqn. @) pro- 
viding the condition needed to fix the relative phase of A and A. 

Let us now consider how the necessary information about the components 
of the amplitude may be extracted from the experimental data. For both the 
uj and p the polarization vector can be completely determined in their rest 
frame by their decays. For p ± — > tc^tt we can define the polarization vector 
as E p = Piro/\P n o\ while for the uj we define E u = P n + x P 7T -/\P 7T + x P n -\. 
Let us introduce 6\ to be the polar angle between E p and the boost axis in 
the p frame and likewise 82 the polar angle between E w and its boost axis. 
In addition we introduce the azimuthal angle ip between E p and E u with 
respect to the p boost axis. It is convenient to describe the final state in a 
helicity bases which contains the following states: \R) = | + l) p | + \L) = 
| — l) w and |0) = |0) p |0) w . In the following we will rewrite the transverse 
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modes as parity eigenstates: \S) = (\R) + \L))/y/2 and \P) = (\R) - \L))/y/2 
where the action of parity is n : |0), \S) -> + |0), \S) and n : \P) -> -|P). 

The amplitudes for production of the |0), (S") and |P) states respectively 
may be expressed in the three component notation: A = [A , As, Ap\. The 
angular distribution (see e.g. [fT^| ) for the two vector decay is dT/d& = 
Y^=i^ift(^) where $ represents the phase space, <i<3> = dzidz 2 d%lj / (87r) and 
the basic distributions are given by: 



h = 9z\z\- j 2 = (9/2)«?«l cos 2 ^; h = (9/2)n 2 M 2 sin 2 ^; 

/ 4 = 9v^2^i^2MiM2 cos^; / 5 = 9v^2-2i sin -0; 

/6 = 9-u 2 m 2 sin ^ cos ip (4) 

with ^ = cos#j and Ui = sin^ (i = 1, 2). In terms of the helicity amplitudes, 
the coefficients are thus: 



X\ — \Aq\ 2 ] X 2 — \A S \ 2 ] X 3 — \A P \ 2 ; , . 

X 4 = Re(A A* s ); X 5 = Im(A A P ); X 6 = Im(A s A* P ); {0> 

we have normalized the coefficients so that Br = X% + X2 + X3. 

A convenient method to determine Xi from the angular distribution is to 
introduce a set of operators [[L3| gi such that < <?j >= Xi/(X\ + X 2 + X 3 ). 
We will use the following set of operators with this property: 



g x = (40/i - 5/2 - 5/ 3 )/126, <? 2 = (265/ 2 - 20/! - 85/ 3 )/504, 

g 3 = (265/ 3 - 20/i - 85/ 2 )/504, ^4,5 = 25/ 4 , 5 /36, ^ 6 = 25/ 6 /72. (6) 

Once an experimental determination of the observables {Xi} has been 
made, it is possible to determine the amplitudes up to the following ambigu- 
ities: (1) An overall phase, (2) the transformation A — > II A*. Thus we have 
six observables determining three complex (amplitudes), or equivalently six 
real parameters. The overall phase ambiguity means that the remaining five 
amplitude parameters are over determined by the six observables. An unob- 
servable phase may be applied to both A and A leaving one only the relative 
phase between the B + and B~ decays undetermined. To be specific, we will 
take this weak phase to be A = a.rg(A (B~) Aq(B + )) . Once A is determined, 
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the experimental observables in B + and B~ decay will allow us to obtain 
bsina from eqn. (|3|). To obtain sin a we therefore need a means to fix A as 
well as information about b. 

Consider now what may be learned from B° — > <fiK*. In the SM we ex- 
pect negligible CP violation in this mode (e.g. in [llj the CP violation is 
estimated to be 0(1%)); thus X 1 _ 4 (0 J fsT* ) = X 1 _ 4 (0^*°) and X 5j6 ((j)K* ) = 
—X 5 fi((j)K* ). Within the SM and using our assumption that the gluon 
does not fragment to a single vector meson, the amplitude A((f)K*°) = 
\f2h(y ts /Vtd) hence the observables will be related to b by: 



\b Q \ 2 = qX x (F), \b s \ 2 = qX 2 (F), \b P \ 2 = qX 3 (F), Re(b b* s ) = qX 4 {F), 
Im(b b* P ) = qX B (f), Im(b s b* P ) = qX 6 (F), (7) 

where JF = <pK* and q = \Vtd/Vt s \ 2 /2. We are now in a position to determine 
A and sin 2 a since eqn. ([F|) gives six equations in only two unknowns (i.e. A 
and a) when combined with eqn. (|||). 

We now estimate the precision with which one may extract sin 2 a using 
this method. The statistical errors will, of course, depend on the number 
of B mesons which are available. To quantify this for a given decay mode 
B — > X, let us define Nb = (number of B)(acceptance for X) for each 
type of 5-meson. For concreteness, we consider Np = 10 8 and 5 x 10 8 in our 
calculations. 

The pure penguin modes B + — > <pK* + and B° — > <f)K*° have in fact 
recently been observed at BaBar @, [0§ with Br(B + -> <pK* + ) = (9.7lg;4± 
1.7) x 10" 6 and Br(B° -> <pK*°) = (8.6+l;| ±1.1) x 1(T 6 . In our sample 
calculation we will take this branching ratio to be 10~ 5 . Estimating the 
corrected b — > d penguin contributions, we get Br penguin (B~ — > p~u)) w 
2\V td /V ts \ 2 Br(B -> <j)K*) w 8x 10" 7 where [|] \V td /V ts \ w 0.2. The tree graph 
for this final state has a similar topology to the observed [20, B + — > 7r + 7r° 
decay modes with branching ratios of ~ 5 ±2 x 10~ 6 . The branching ratio to 
a two vector final state should be somewhat larger because of the additional 
helicity states. Here we will assume that this will increase the branching ratio 
by roughly a factor of 3 as is the case for Br(B + — > p + D*°) / Br(B + — > tt + D°) 
so that Br tree (B + -> p+uu) w 15 x 1(T 6 @. 

The precision for extracting sin 2 a will depend on the distribution of the 
amplitudes between the various helicity states. As an example, which we will 
call Case 1, we will take the toy model for the amplitudes given by: 



7 



a oc [1,1,-1] b oc [0.23,0.23,0.23] 



(8) 



which has the correct ratio of total rates. In this toy model the amplitudes 
have no relative phase from final state interactions. Case 2 will be the corre- 
sponding results averaged over helicity distributions and phases subject only 
to the condition that |a| and |b| are fixed. 

To simulate the reconstruction of sin 2 a from experimental data we will 
calculate the observables for each mode and estimate the correlated errors 
assuming that these are measured using the operators </j. The tree-penguin 
mode give us A and A up to a relative phase. We can thus use eqn. to 
get b sin a so that eqn. fl7|) further constrains the fit and a minimum value of 
X 2 can be obtained for each value of sin 2 a. The dependence of \ 2 on s i n2 a 
thus allows us to obtain the statistical error A sin 2 a given in Table |2] where 
we take a = 30° and 60°. In the last line of this Table, we consider a similar 
exercise for a mode sensitive to 7. 

Table ||] gives results for the tree penguin modes: B + — > p + u where 
the tree is color allowed, B° — > p°u where the tree is color suppressed and 
a combined fit with both modes taken together. We also consider the B s 
mode, B s —> K* + p~. In each case we take Nb = 10 8 and 5 x 10 8 and give 
results for Cases 1 and 2. 

Let us now consider how the validity of the approximations used in this 
method may be experimentally verified. The validity of SU(3) for penguins 
may be tested by comparing the various pure penguin modes considered. If 
this symmetry is exact, the magnitude and relative phases of all the b — > s 
pure penguin modes should be the same. 

The assumption that EWPs do not play a large role is required in most 
isospin analysis ||. [7|. § since the EWPs have the isospin properties of the tree 
but the weak phase of the penguin so it cannot be cleanly separated from the 
tree without additional information. Its presence may be specifically checked 
in the decay B^ — > p ± p°. The final state must have total isospin / = 2 so in 
the SM only the tree and EWP can contribute. A significant EWP contri- 
bution which has a relative weak phase of a with respect to the tree would 
thus lead to CP violation. There are six possible CP violating observables 
that can be obtained from the distributions: Xi__ 4 (p + p°) — Xi_ 4 (p _ p°) and 
X5fi(p + p°) + X 5i6 (p~p ). Note that the first four are even under naive time 
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Table 2: The error in sin 2 a or sin 2 7 assuming the central value of a, 7 = 
30° or 60°. The results under Case 1 use the helicity distribution given in 
eqn. (|8|) while Case 2 is a Monte Carlo average both over helicity distribution 
and over phases keeping the branching ratios fixed. In each entry the first 
number is for Ng = 10 8 while the second is for Nb = 5 x 10 8 . 



Tree- Penguin 


CKM quantity 


Cat 

30° 


;e 1 

60° 


Cas 

30° 


,e 2 

60° 


B + -> p+LU 


sin 2 a 


.065 


.029 


.112 


.050 


.063;. 033 


.115;.062 


B° -> p u 


sin 2 a 


.220 


.118 


.221 


.107 


.095; .051 


.114;.060 


B + -> p + tu & B° -> 


sin 2 a 


.056 


.025 


.101 


.045 


.055;. 023 


.081;.040 


fi s -► 


sin 2 a 


.088 


.039 


.151 


.068 


.098;. 039 


.146;.072 


5+ -> K* + p- 


sin 2 7 


.054 


.032 


.069 


.049 


.099; .054 


.095; .057 



reversal fl23j, T N , thus require a rescattering phase while the last two are T N 



odd and do not. The analogous test for the EWP can also be performed in 
the decay B ± — > tt^tt where an EWP could generate a partial rate asym- 
metry. The possible advantage of using the p ± p° final state is that there are 
six CP odd signals rather than one in 7r ± 7r°. 

In the Kn, K*tt, Kp and K*p systems one can also check for EWP 
effects by looking for CP violation with isospin properties inconsistent with 
tree penguin interference 



14, 24 1 . In the Kti system the absence of EWP 



implies: 



2A(£- -> K-u Q ) - A(B~ -> K\-) - A(B -> K-ji + ) + 2A(B -> K°p°) = (9) 

where A(X) is the partial rate asymmetry for the reaction X. Eqn. (|9]) also 
applies to Kp and K*tt while for K*p it applies separately to X±, X2 and 
X 3 . 

It is also possible to reduce the potential impact of the EWP by choosing 
modes in which their contribution will be color suppressed. In particular, the 
tree penguin interference processes B s — > K* + p~ and B — > p + p~ will color 
suppress the EWP as will the pure penguin process B s —>■ K*°K*°. 
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The assumption that the gluon does not fragment into a single vector 
meson is true at lowest order in perturbation theory by color conservation. 
We can estimate how well it holds when QCD corrections are introduced, 
through renormalization group improved perturbation theory. If we define 
b' to be the penguin contribution where the quarks from the gluon form a 
single meson, the ratios of the contributions should be: 



|b'|/|b| « (3C 3 + C 4 + 3C 5 + C 6 )/(C 3 + C 4 + C 5 + C 6 ) (10) 



where we have used the notation of [25) . With the numerical results from that 



paper we find that this ratio is about 0.04 for various values of Aqcd both 
in naive dimensional regularization and in the 't-Hooft — Veltman schemes 
at next to leading order. 

To summarize, we have shown that unitarity angles a and 7 can be ex- 
tracted by a method that is rather unique in that it only uses direct CP 
violation and does not need any time dependent measurements. We gave 
several examples of modes which can be used and also discussed how ex- 
perimentally the underlying theoretical approximations can be tested. Since 
these angles are fundamental parameters of the SM, and their precise deter- 
minations could lead to its crucial tests, the importance of measuring them 
in several different ways can hardly be over emphasized. 

This research was supported in part by US DOE Contract Nos. DE-FG02- 
94ER40817 (ISU); DE-AC02-98CH10886 (BNL). 
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